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ABSTRACT 

We have investigated the properties of a group of 7-ray emitting globular 
clusters (GCs) which have recently been uncovered in our Galaxy. By correlating 
the observed 7-ray luminosities with various cluster properties, we probe the 
origin of the high energy photons from these GCs. We report is positively 
correlated with the encounter rate Tc and the metalicity [Fe/H] which place 
an intimate link between the gamma-ray emission and the millisecond pulsar 
population. We also find a tendency that L-y increase with the energy densities 
of the soft photon at the cluster location. Furthermore, the two-dimensional 
regression analysis suggests that L^, soft photon densities, and rc/[Fe/H] possibly 
span fundamental planes which potentially provide better predictions for the 7- 
ray properties of GCs. 

Subject headings: gamma rays: stars — globular clusters: general — pulsars: 
general 



1. INTRODUCTION 

Millisecond pulsars (MSPs) are generally believed as the descenders of the low-mass 
X-ray binaries (LMXBs) (Alpar et al. 1982). As the formation rate per unit mass of LMXBs 
is orders of magnitude greater in globular clusters (GCs) than in the Galactic field (Katz 
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1975; Clark 1975), it is not surprise that 80% of the detected MSPs are located in Gcfl (cf. 
Manchester et al. 2005). The relatively high formation rate of LMXBs and MSPs is a natural 
consequence of the frequent stellar encounters. With the X-ray populations in various GCs 
have been revealed by the Chandra X-Ray Observatory, Pooley et al. (2003) and Gendre et 
al. (2003) have found a positive correlation between the number of LMXBs in GCs and the 
stellar encounter rate, Fc. This provides evidence for the dynamical formation of LMXBs 
in GCs. As the descenders of the LMXBs, MSPs are also expected to have a dynamically 
origin. 

Very recently, with the corrections of the observational effects in the radio pulsar surveys 
toward different GCs, Hui, Cheng & Taam (2010) have found a positive correlation between 
the MSP populations in GCs and T^, which has long been predicted. Moreover, the authors 
have also found another positive correlation between the metalicity and the MSP population. 
This relation is not unexpected as the high metalicity in a GC can result in a more efficient 
orbital shrinkage by magnetic braking. Therefore, the parameter space for the successful 
Roche-lobe overflow is enhanced (Ivanova 2006) and subsequently lead to a higher formation 
rate of MSPs. 

A brand new window for investigating the MSPs in GCs has been open by the launch 
of the Fermi Gamma-ray Space Telescope. Since MSPs are the only steady 7-ray emitters 
in GCs, observations with Fermi can provide an alternative channel for investigating MSP 
populations. Shortly after the commence of its operation, the Large Area Telescope (LAT) 
onboard the spacecraft has detected 7— rays from 47 Tucanae (hereafter 47 Tuc) (Abdo et al. 

2009) . Terzan 5, which hosts the largest known MSP population, has also been subsequently 
detected (Kong et al. 2010). As the sensitivity of LAT increases monotonically with the 
continuous all-sky survey, a total of 15 confirmed detections of 7— ray emitting GCs have 
very recently been reported by Abdo et al. (2010a; 2010b) and Tam et al. (2010). Using 
~ 1.5 years of LAT data, Abdo et al. (2010a) have detected 6 new 7-ray GCs besides 47 Tuc 
and Terzan 5. On the other hand, Tam et al. (2010) have recently reported 7 other new 
detections with ~ 2 years data. Among all known 7-ray GCs, Liller 1, which has the highest 
metalicity in our Milky way, is also found to have the highest 7-ray luminosity (Tam et al. 

2010) . This discovery further suggests that the effects of metalicity cannot be neglected. 
Thanks to these surveys, we are able to study these clusters as a unique class for the first 
time. 

To explain the 7— rays from GCs, there are two main streams. Venter & de Jager (2008) 
and Venter et al. (2009) suggest the 7— rays are originated from the curvature radiation of 
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electrons in MSP magnetospheres. However, the frequent stellar interactions can lead to a 
complicated magnetic field structure (Cheng & Taam 2003), which can possibly explain the 
difference of the radio and X-ray properties of MSPs in GCs with respect to those located in 
the Galactic field (cf. Hui, Cheng & Taam 2009, 2010). One consequence of the complicated 
surface magnetic field is to turn off the accelerating region for producing high energy photons. 
Ruderman & Cheng (1988) argue that if the open field lines are curving upward due to the 
effect of local field then in this case e~ /e^ pair production and outflow can occur on all 
open-field lines. As a result, the outer- magnetospheric gap is quenched by these pairs. This 
scenario is supported by the fact that the MSPs in 47 Tuc are essentially thermal X-ray 
emitters (Bogdanov et al. 2006). All these demonstrate the potential difficulties of the 
pulsar magnetospheric model in explaining the observed 7— rays from GCs, which motivate 
the exploration of additional / alternative emission mechanisms. 

On the other hand, Bednarek & Sitarek (2007) have proposed that the inverse Compton 
scattering (ICS) could be a possible mechanism to produce 7— rays from GCs. In their 
model they predict that GCs could be sources of GeV-TeV photons. Unfortunately they 
have ignored the contribution from the Galactic background photons. Very recently, through 
generalizing the ICS model by including various soft photon fields, Cheng et al. (2010) have 
found that the observed 7— ray spectra of all 8 7-ray GCs reported by Abdo et al. (2010a) 
can be well- modeled by the ICS between relativistic electrons/positrons in the pulsar wind of 
MSPs in the GCs and the background Galactic soft photons. This provides another possible 
explanation for the origin of the 7— rays. 

In this paper, we report the results from exploring the 7— ray emission properties by 
comparing with various cluster properties, which provides us with insight on the origin of 
the 7— rays from this class of GCs. In §2, we report the method and the results from the 
correlation and the regression analysis. We subsequently discuss the implication of these 
results in 53. 



2. CORRELATION & REGRESSION ANALYSIS 

Abdo et al. (2010a) and Tam et al. (2010) have reported 15 GCs with firm detec- 
tions. This sample size is somewhat larger than that adopted by Pooley et al. (2003) for 
investigating the relations between the X-ray point source populations in GCs and various 
cluster parameters. The sensitivity limit of the current sample is ~ 6 x 10^^^ erg cm~^ s~^ 
(0.1-100 GeV). The properties of 15 confirmed 7— ray GCs are summarized in Tableland 
the entries are explained in the following. 
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Table 1. Properties of the 7— ray emitting GCs. 



Cluster Name 


kpc 


r ^ 

c 


[Fe/H]^ 




^optical 

eV cm~^ 


eV cm~^ 


log L^^ 
erg s~^ 


Adopted from Abdo ct al. (2010a) 


47 Tuc 


4.0 


44.13 


-0.76 


-9.17 


0.93 


0.25 


34.681}^ 


Omega Cen 


4.8 


4.03 


-1.62 


-10.07 


1.61 


0.51 




M 62 


6.6 


47.15 


-1.29 


-9.09 


8.07 


0.86 




NGC 6388 


11.6 


101.99 


-0.60 


-9.74 


2.59 


0.56 


35.4ll°i^ 


Terzan 5 


5.5 


118.29 


0.00 


-6.51 


7.02 


1.37 


35AltrZ 


NGC 6440 


8.5 


74.17 


-0.34 


-8.78 


10.79 


1.00 


35.381°-}^ 


M 28 


5.1 


13.10 


-1.45 


-7.98 


5.47 


0.92 


3A.79t'oil 


NGC 6652 


9.0 


1.24 


-0.96 


-6.43 


3.65 


0.51 


34.89l°:i« 


Adopted from Tam et al. (2010) 


Liller 1 


9.6 


77.98 


0.22 


-7.63 


10.53 


1.40 


77+0.13 
"J^- ' ' -0.18 


M 80 


10.3 


31.31 


-1.75 


-8.23 


1.88 


0.33 


34.92+°-28 


NGC 6441 


11.7 


88.42 


-0.53 


-9.64 


3.59 


0.69 


35.571°:?^ 


NGC 6624 


7.9 


14.65 


-0.44 


-7.49 


6.03 


0.67 


35.17+°?? 


NGC 6541 


6.9 


20.00 


-1.83 


-8.34 


4.69 


0.61 


34.541°!^ 


NGC 6752 


4.4 


10.78 


-1.56 


-7.94 


2.01 


0.48 


34.141^:^ 


NGC 6139 


10.1 


13.28 


-1.68 


-8.36 


4.10 


0.69 


35.03_o:34 



'^Cluster distance adopted from Abdo et al. (2010a) and Tam et al. (2010). 

^Two-body encounter rate estimated by Por^a^^ with the value scaled with that in M4 which 
has po = 10^-^^ Lqpc~^, Tc = 0.53 pc and ao = 8.9 km/s 

^Metalicity 

^Absolute visual magnitude 

^Energy densities of various soft photon fields (see text) 

7— ray luminosities adopted from in Abdo et al. (2010a) and Tam et al. (2010) 
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For choosing the cluster parameters for the correlation analysis, we have considered 
two-body encounter rate Fc, metalicity [Fe/H], absolute visual magnitude My, as well as 
Galactic background optical / infrared photon densities at the locations of the GCs ^optical 

/ «IR- 

Fc is the most obvious parameter related to the binary formation rate and hence the 
number of MSP in a GC. This parameter can be estimated as P^t^cTq^ where po is the 
central luminosity density, is the core radius and do is the velocity dispersion at the 
cluster center, ctq are adopted from Gnedin et al. (2002). For po and rc, the values are taken 
from Harris (1996; 2003 version) and modified for the distances adopted for this analysis 
(cf. Tab. [1]). Besides Fc, Hui et al. (2010) have shown that [Fe/H] are also a key parameter 
in determining the size of the MSP population in a GC. The values of [Fe/H] are taken 
from Harris (1996). On the other hand, if stellar encounters were not the major channel 
of the binary formation, one would expect the binary population to be correlated with the 
cluster mass Mqq for a primordial binary origin (Lu et al. 2009; Lan et al. 2010). Pooley 
et al. (2003) have estimated the cluster mass by integrating the King's profiles of the GCs. 
And therefore, their mass estimates are naturally correlated with F^. On the other hand, 
assuming a constant mass-to-hght ratio, Mqq can also be estimated from the absolute visual 
magnitude My: Mqc = lO""-^^^ (cf. Hui et al. 2010; Lu et al. 2009). Different from the 
estimates adopted by Pooley et al. (2003), the correlations between our mass estimates with 
Fc is only confident at the level less than 53% (see Figure [1^). Also, the correlation between 
My and [Fe/H] only attains a confidence level < 56%. The values of My are also taken from 
Harris (1996) and modified for the adopted distances as presented in Tabled! 

Apart from the number of MSPs, in the context of ICS model, the 7— ray luminosity of 
a GC also depends on the energy density of the soft photon field (see Cheng et al. 2010). 
There are three components of background photons in the Galaxy which can interact with 
the relativistic leptons: they are relic, infrared and optical photons. As the energy density 
of the relic photons is uniform and does not vary from cluster to cluster, we ignore it in our 
analysis. We obtain the estimates of Galactic optical and infrared photon density, ^optical 
and uiR with the GALPROP code (Strong & Moskalenko 1998). 

Without a priori knowledge of the distributions of the tested quantities, we follow Pooley 
et al. (2003) to adopt a nonparametric correlation analysis. The computed Spearman rank 
correlation coefficients between and various tested quantities are tabulated in Table El 
We have also computed the linear correlation coefficients (i.e. Pearson's r) for an intutitive 
account for the data scattering, though they are less robust than the Spearman ranks in 
quantifying the correlations. We have also performed the 1-dimensional regression analysis. 
The best-fit parameters are also given in Table [2j All the quoted uncertainties are 95% 
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confidence intervals. The best-fit relations of these quantities with are plotted as solid 
straight lines in Figure [T|d and Figure [2j We have also plotted the upper and lower 95% 
confidence bands for a visual comparison for the data scattering in each panel. 

Among all these parameters, the weakest correlation is found for the log — My rela- 
tion which is only significant at ~ 15% confidence level. Therefore, there is no convincing 
correlation between these two quantities (see Figure [T|d). On the other hand, the correlations 
of with Fc and [Fe/H] are confident at a level over 99% (see Figure H]). All these findings 
are fully consistent with the results from analysing the radio MSP population (Hui et al. 
2010). 

While the correlation between and Fc was reported by Abdo et al. (2010a) with 
8 GCs, the effect of metalicity was ignored in their work. From the 95% confidence bands 
shown in Figure |21 the degree of data scattering of the logL^ — [Fe/H] relations is found to 
be the smallest among all the tested single parameters. This can be also reflected by the 
fact that its corresponding linear correlation attains the confidence levels over 99.9%. 

For the tested soft photon fields. Figure |2] shows that also tends to increase with their 
energy densities. For Wopticai and uir, the confidence levels for the correlations are > 96% 
and > 99% respectively. 

In comparison to the logL^ — logF^ and logL^ — [Fe/H] relations, the relatively large 
scattering of the data points in the plots of log — log Wopticai and log — log mjr can be due 
to the fact that the GALPROP code is a simplification of the real situation in the Galaxy, 
which assumes an axisymmetric distribution of all parameters of the program (cf. Strong & 
Moskalenko 1998). Although it provides reliable average density of background photons in 
the Galaxy, the real photon densities at the locations of the GCs can differ by a factor of a 
few, in particular for those close to the disk. 

Despite the scattering, the correlation analysis strongly suggests is likely related to 
the soft photon energy density estimates. This inference is consistent with the ICS model- 
ing the 7— ray spectra of GCs (Cheng et al. 2010) which indicates neither the number of 
MSPs nor the soft photon energy density is the sole factor in determining L^. With this 
consideration, we investigate if L^, rc/[Fe/H], and Wopticai/wiR span a fundamental plane 
by a 2-dimensional regression analysis. We have examined the sample with the following 
relations: 



log L^ = ai + 02 log Fc + as log Wopticai 



(1) 



7 



log L. 



7 



04 + as log + fle log uir 



(2) 



log = ay + ag [Fe/H] + ag logWoptkai 



(3) 



log = aio + an [Fe/H] + au logWiR. 



(4) 



The best-fit parameters are tabulated in Table [31 We have shown the edge-on view of 
these best-fit fundamental plane relations in Figure [31 In comparison with Figure [H the 
data scatter in these plots are somewhat reduced which suggest these fundmental plane 
relations can possibly provide us with better 7— ray luminosity predictors than the single 
parameter relations. To better constrain the uncertainties of these parameters, we have 
further computed the la, 2a and 3a confidence contours for various parametric spaces which 
are shown in Figure [H 



We have examined the 7— ray emission properties of a group of GCs. By investigating 
the possible correlations between the 7— ray power and a number of cluster properties, we 
shed light on the origin of the 7— rays from these GCs. First of all, the correlation between 
and suggests the high energy radiation are intimately related to the population of 
dynamically-formed objects, which are presumably MSPs, confirming Abdo et al. (2010a) 
who used 8 GCs in their study. Together with the lack of any correlation with My and hence 
the cluster mass, this is fully consistent with the inference suggested by Hui et al. (2010) 
and consolidates the dynamical formation scenario of MSPs in GCs. 

Apart from F^, we have found that is also positively correlated with [Fe/H]. This is 
well-consistent with the tendency deduced from studying the radio MSP population in GCs 
(Hui et al. 2010) and the fact that the GC possesses the highest [Fe/H] also has the highest 
(Tam et al. 2010). Ivanova (2006) proposes that the absence of the outer convective 
zone in metal-poor main sequence donor stars in the mass range of O.SSM© - 1.25M0, in 
comparison to their metal rich counterparts can be responsible, since the absence of magnetic 
braking in such stars precludes orbital shrinkage, thereby, significantly reducing the binary 
parameter space for the production of bright LMXBs. For a conventional scenario, MSPs 
are the old pulsars that have passed through the death-line in P — P diagram which are 



3. 



DISCUSSION 
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Table 2. Correlation and 1-D regression analysis of log versus various cluster 

properties. 



Parameters 


Spearman rank 


Probs^ 


Pearson's r 


Probp'^ 


rrf 




logr. 


0.7918 


0.9996 


0.6414 


0.9900 


0.50 ±0.16 


34.38 ± 0.25 


[Fe/H] 


0.7614 


0.9990 


0.7912 


0.9996 


0.59 ±0.15 


35.56 ±0.15 


Mv 


-0.0536 


0.1505 


-0.0767 


0.2141 


0.04 ± 0.08 


35.40 ± 0.67 


log ^optical 


0.5523 


0.9672 


0.5976 


0.9814 


0.78 ±0.27 


34.62 ±0.19 


log um 


0.6589 


0.9925 


0.5970 


0.9812 


1.29 ± 0.44 


35.33 ±0.12 



^The probability that the Spearman rank correlation coefficient is different from zero. 

^The probabihty that the hnear correlation coefficient (i.e. Pearson's r) is different from 
zero. 

'^The best-fits for log — mx ± c where x is the corresponding parameters listed in 
column 1. 



Table 3. Best-fit fundamental plane relations of 7— ray GCs. 



Parameters 


Best-fit values 




34.12 ±0.29 


02 


0.42 ±0.17 


03 


0.62 ±0.29 


04 


34.70 ± 0.30 


05 


0.39 ±0.18 


oe 


0.96 ± 0.49 


a^ 


35.21 ±0.29 


as 


0.49 ±0.18 


ag 


0.44 ±0.31 




35.61 ±0.16 


an 


0.48 ±0.17 


ai2 


0.76 ±0.50 
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subsequently spun-up in the binaries. As the metahcity determines the parameter space for 
successful Roche-lobe overflow, it is also a key parameter in determining the intrinsic number 
of MSPs in a GC (Hui et al. 2010; Ivanova 2006). 

We note that the link between the LMXBs in extragalactic GCs and the metalicity is 
somewhat weaker than with the cluster mass (e.g. Sivakoff et al. 2007; Kim et al. 2006; 
Kundu et al. 2002), which is different from the inference drawn from our investigation of 

the Galactic MSP-hosting or 7— ray selected clusters. However, a direct comparison between 
these two populations has to be cautious. As the MSPs are long-lived and are produced by 
the previous generations of LMXBs, their dynamical properties might be different from that 
of the LMXB population currently observed. Since the relaxation time at the cluster core is 
generally longer than the hfetime of LMXBs, the cluster is continuously evolved with mass 
segregation at the cluster center which can result in a varying formation rate of compact 
binaries (cf. the discussion in Hui et al. 2010). Also, while a large number of LMXB-hosting 
GCs in Virgo cluster early-type galaxies have relaxation times > 2.5 Gyr, there is no single 
GC in our Galaxy with a relaxation timescale larger than this value contains an active LMXB 
(cf. Sivakoff et al. 2007). Although the reason is still unclear, this suggests possible different 
properties between the Milky Way GCs and the extragalactic ones. Further investigations 
are required to understand the difference. 

It is instructive to compare the fundamental plane relations of the 7— ray population 
with the best-fits inferred from the radio MSP population in GCs. Hui et al. (2010) have 
found that the slopes of logFc and [Fc/H] inferred from the radio GC MSPs population are 
0.69 ± 0.11 and 0.72 ± 0.11 respectively Within these quoted uncertainties, the slope of 
the [Fe/H] relation for the radio population is found to intersect with the 2(7 error contours 
for the corresponding parameters of the 7— ray fundamental plane relations (i.e. Og and 
Oil). On the other hand, the logarithmic slope of the Fc relation for the radio population 
is only marginally overlapped with the rims of the 3cr error contours for the corresponding 
parameters inferred from the 7— ray population (i.e. 02 and 05). 

We have also identified possible positive correlations with various soft photon fields 
which have significances compatible with those for the encounter rate and the metalicity. 
These correlations are not expected from the magnctosphcric model. Together with the 
uncertainty of the sustainability of the outergaps in the MSPs in GCs (see §1), our finding 
further motivate the exploration of alternative explanations for the origin of the observed 
7-ray from GCs. 

Abdo et al. (2010a) have argued that the 7-ray emission is magnetospheric in nature 
because of the hard photon indices and the cutoff energies inferred from the phenomenolog- 
ical model is consistent with the values expected from the magnetospheric model. On the 
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other hand, Cheng et al. (2010) have recently found that the ICS model can also describe the 
observed 7— ray spectra of all the GCs discovered by Abdo et al. (2010a) very well. Simply 
based on the model fitting, we were not able to discriminate these two scenarios unambigu- 
ously. However, different from the case of the magnetospheric model, positive correlation 
between the energy density of the soft photon fields are expected in a ICS scenario as the 
ICS power is directly proportional to soft photon energy density. 

The energy density of the background soft photon field depends on the location of the 
cluster. We notice these 7-ray GCs arc possibly resided in the Galactic bulge and therefore 
they arc also mctal-rich clusters. This results in a natural correlation between the mctalicity 
and the soft photon energy density with a significance > 95%. Therefore, it is non-trivial to 
disentangle the effects of these two parameters. 

In any case, our investigation strongly suggests that either the metalicity or the soft 
photon energy density has to be the new parameter, in addition to Fc, in determining the 
observed 7-ray luminosities. This inference is supported by comparing the results reported 
by Abdo et al. (2010a) and Tam et al. (2010). Apart from the 8 confirmed cases, Abdo 
et al. (2010a) have also reported 5 non-detections which include three upper-limits and two 
other cases with the 7— ray emission shghtly offset from the respective GC cores. With the 
LAT data of a longer exposure, Tam et al. (2010) have found a larger number of 7-ray GCs 
including 4 previously non-detected cases in Abdo et al. (2010a). This leaves M 15 to be 
the only non-detected GC in Abdo et al. (2010a). This is not unexpected from trends of 
mctalicity and background soft photon energy density. Although the encounter rate of M 15 
(Fc = 53.9) is even higher than M 62, its metalicity ([Fe/H]=-2.26) and the background soft 
energy densities at its location (^optical ~ 0.44 eV cm"^; mir = 0.11 eV cm~^) are much lower 
than those of all 15 confirmed 7-ray GCs. 

In view of the aforementioned complication, it is not possible to discriminate the ICS 
and the magnetospheric scenarios unambiguously simply based on the currently available 
information. Also, there is still a degeneracy within the context of ICS model. Cheng et 
al. (2010) show that the 7— ray spectrum from 47 Tuc can be explained equally well by 
upward scattering of either the relic photons, the Galactic infrared photons or the Galactic 
optical photons whereas the 7— ray spectra from the other seven GCs reported by Abdo et 
al. (2010a) are best fitted by the upward scattering of either the Galactic infrared photons 
or the Galactic optical photons. This has prompted us also to discriminate which source 
provides the predominant soft photon field for ICS. 

Since the IC radiation power is directly proportional to the energy density of the soft 
photon field, a logarithmic slope of unity is thus expected for the fundamental plane param- 
eters 03, oe, Qg and 012. For the corresponding parameters of mir (i.e. a% and 012), the line of 
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unity is found to cut through the centers of the la error contours for both Fc and [Fe/H] fun- 
damental plane relations. On the other hand, for the parameters correspond to i^opticai (i-e. 
03 and ag), the expected value is only marginally intersected with their 3a error contours. 
Although an unambiguous conclusion cannot be drawn from the current population yet, the 
comparison between the theoretical expectation and the fundamental plane parameters docs 
favor the scenario involving the background infrared emission as the soft photon field for IC 
upscattering. 

We would like to point out that the predicted spectral shape in energy regime much 
larger than 10 GcV is significantly different for different soft photon fields (cf. Fig. 2 in 
Cheng et al. 2010). Therefore, observations with TeV facilities, such as MAGIC, HESS and 
VERITAS, can be feasible to lift up this degeneracy. Furthermore, as the ICS model and 
the magnetospheric model predict a rather different TeV spectrum for GCs (Cheng et al. 
2010; Venter et al. 2009), TeV observations in the future can possibly better discriminate 
these two possible contributions of soft photons. 

Constraints for the emission model can also be derived from the other energy bands. 
With the diffusion of the relativistic pulsar wind particles, it has been shown that extended 
radio and X-ray emission from the GCs can also be produced by synchrotron radiation and 
ICS respectively (Cheng et al. 2010). This is consistent with the recent discovery of the 
diffuse X-rays around Terzan 5 which exteneded up ~ 10 pc (Egcr, Domainko & Clapson 
2010). Although a clear scenario cannot be identified yet, these diffuse X-rays arc more likely 
to have a non-thermal origin (Egcr, Domainko & Clapson 2010). Assuming these X-rays are 
originated from the tail of ICS, the corresponding 7— ray spectrum can be calculated (Cheng 
et al. 2010). Therefore, a systematic search for the extended X-ray and radio feature outside 
the half-mass radii of the other 7— ray GCs can provide us indpendent constraints. 

In exploring the fundamental plane relations, our analysis suggests that by combining 
the soft photon energy densities with Fc/fFe/H] the data scattering can be reduced. These 
best-fit relations can provide the indicators in identifying what kind of GCs are potential 
7— ray sources for a further search. And the other way round, any deeper 7— ray search from 
the GCs can result in an enlarged sample size and a lower sensitivity limit than the current 
value (i.e. 6 x 10~^^ erg cm~^ s~^), which will certainly enable a further test for all these 
reported relations. 
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Fig. 1. — a. Two-body encounter rate Tc vs. absolute visual magnitude My. b. The 7— ray 
luminosity vs. My- The straight line represents the best-fit straight line with the errors 
of the data points fully taken into account. The dotted lines represent the upper and the 
lower 95% confidence bands. 
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Fig. 2. — vs. various individual cluster properties. The straight hues in the plots 
represent the best-fits from the linear regression with the errors of the data points fully 
taken into account. The dotted lines represent the upper and the lower 95% confidence 
bands. 
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Fig. 3. — The edge-on views of the fundamental plane relations of 7— ray GCs. The straight 
lines in the plots represent the projected best-fits given in Table El 
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Fig. 4. — The maps for various parametric spaces of the fundamental plane relations. The 
dashed lines illustrate the la, 2a and 3a confidence contours for two parameters of interest 
which encircle the best-fit values (i.e. the positions with the lowest x^)- 



